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ABSTRACT 

We study the properties of Brightest Cluster Galaxies (BCGs) drawn from a catalogue 
of more than 69000 clusters in the SDSS DR6 based on the adaptive matched filter 
technique (AMF, Szabo et al., 2010). Our sample consists of more than 14300 galaxies 
in the rcdshift range 0.1-0.3. We test the catalog by showing that it includes well-known 
BCGs which lie in the SDSS footprint. We characterize the BCGs in terms of r-band 
luminosities and optical colours as well as their trends with redshift. In particular, 
we define and study the fraction of blue BCGs, namely those that are likely to be 
missed by either colour-based cluster surveys and catalogues, as shown by a direct 
comparison to maxBCG clusters that are matched in the Szabo et al. catalogue. We 
further compare the properties of the BCGs to those of the second and third brightest 
galaxies in the same cluster. Finally, we morphologically classify those galaxies hosted 
in the richest clusters. 

We find that the BCG luminosity distribution is close to a Gaussian, whose mean 
has a redshift evolution broadly consistent with pure aging of the galaxies. Richer 
clusters tend to have brighter BCGs, however less dominant than in poorer systems. 
4-9% of our BCGs are at least 0.3 mag bluer in the g-r colour than the red-sequence 
at their given redshift. Such a fraction decreases to 1-6% for clusters above a richness 
of 50, where 3% of the BCGs are 0.5 mag below the red-sequence. A preliminary mor- 
phological study suggests that the increase in the blue fraction at lower richnesses may 
have a non-negligible contribution from spiral galaxies. In terms of redshift evolution, 
the overall blue fraction goes from ~5% in the redshift range 0.1-0.2 to ~10% in the 
rcdshift bin 0.2-0.3. The blue fraction seems to increase at higher redshifts, however 
the scatter in the colours and the fact that the catalog is no longer complete hamper 
us from having firm conclusions. We show that a colour selection based on the g-r 
red-sequence or on a cut at colour u-r> 2.2 can lead to missing the majority of such 
blue BCGs. Finally, the blue fraction increase by a factor 1.5 at most when the study 
is extended to the three brightest galaxies of each cluster. 

We also extend the colour analysis to the UV range by cross-matching our cata- 
logue with publicly available data from Galex GR4 and GR5. We show a clear corre- 
lation between offset from the optical red-sequence and the amount of UV-excess. 

Finally, we cross-matched our catalogue with the ACCEPT cluster sample (Cav- 
agnolo et al., 2009), and find that blue BCGs tend to be in clusters with low entropy 
and short cooling times. That is, the blue light is presumably due to recent star for- 
mation associated to gas feeding by cooling flows. 

Key words: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - 
galaxies: evolution - cooling flows - X-rays: galaxies: clusters 



1 INTRODUCTION 
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The Brightest Cluster Galaxies (BCGs) are the most mas- 
sive galaxies in the universe, with most of their stellar mass 
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in place by redshift 2. Therefore, they are expected to experi- 
ence the galaxy formation process in the most extreme way. 
Namely they should form at earlier times, more rapidly and 
with a more intense star formation event than lower mass 
cluster members. While such a formation scenario is natu- 
rally explained in the framework of the revised monolithic 
collapse models (e.g. Larson, 1974, Pipino et al., 2008), it 
seems more difficult to reconcile their existence within the 
hierarchical growth scenario where the largest structures are 
the last to form. A possible way out is that feedback halted 
the star formation at very early times and the mass assembly 
of BCGs can be simply explained with a series of gas-less 
mergers of old stellar systems (De Lucia & Blaizot, 2007, 
but see Whiley et al., 2008, Pipino et al., 2009b, Pipino & 
Matteucci, 2008). A closer look at their stellar populations 
tells us that BCGs have similar mean stellar ages and metal- 
licities to non-BCGs ellipticals of the same mass but they 
have somewhat higher a/Fe ratios, indicating that star for- 
mation may have occurred over a shorter time-scale in the 
BCGs (von der Linden et al., 2007, Loubser et al. 2008). 
Moreover, they depart from the Faber-Jackson relation for 
ellipticals (Faber & Jackson, 1976) and seem to be larger 
and with a higher stellar velocity dispersion than ellipticals 
of the same mass (von der Linden et al., 2007, Bernardi et 
al., 2007, 2008). Furthermore, the BCG luminosity function 
differs from the usual Schechter (1976) form that holds for 
normal cluster members (e.g. Hansen et al., 2005), in that it 
can be modelled as a Gaussian whose mean increases with 
the cluster richness (Lin et al., 2004). Therefore, a study of 
the colours in BCGs as opposed to "ordinary" early-type 
galaxies out to high redshift is a test bench to discriminate 
among models for galaxy formation (e.g. Roche et al., 2009). 
In particular, the colour evolution may place constraints 
to the time and the intensity of their last star formation 
episode, whereas a joint analysis of the colour evolution with 
the thermal status of the surrounding intracluster medium 
may tell us what halted the gas supply and inhibited further 
star formation. 

Luminous red galaxies - as most of the BCGs are - are 
used as a good tracer of the underlying dark matter distri- 
bution (Ho et al., 2009a, Reid & Spergel, 2009). Their prop- 
erties are strongly linked to the host halo mass, so their cen- 
sus can provide us with a map of large-scale over-densities 
in the Universe. They have been used to detect baryonic 
acoustic oscillations (e.g. Sanchez et al, 2009, and references 
therein), and it has been put forward (Ho et al., 2009b) that 
a combination of a galaxy redshift survey such as SDSS and 
a CMB survey can be used as a method for detecting the 
missing baryons. Therefore BCGs can be a promising tool 
for precision cosmology as well. 

Finally, since BCGs occupy a special place in that they 
sit at the bottom of the cluster potential well, we expect 
their present-day properties to be linked to the state of the 
intracluster gas. Recent studies have reported examples of 
ongoing star formation in the massive central galaxies of cool 
core clusters (Crawford et al. 1999, Edge 2001, Goto 2005, 
McNamara et al. 2006, Hicks & Mushotzky 2005, O'Dea et 
al., 2006, Edwards et al., 2007, 2009). Bildfell et al. (2008) 
found that the presence of optical blue cores in 25% of its 
BCG sample is directly linked to the X-ray excess of the host 
clusters. Moreover the position of these BCGs coincide with 
the peak in X-ray emission. Their interpretation is that the 



recent star formation in BCGs is associated with the bal- 
ance between heating and cooling in the ICM in the sense 
that the clusters that are actively cooling are forming stars 
in their BCGs. Other evidence comes from Rafferty et al 
(2008), Cavagnolo et al. (2008). In particular, Pipino et al. 
(2009a) demonstrated a one-to-one correspondence between 
blue cores in BCGs and a UV-enhancement observed using 
GALEX. The blue light coming from the cores might ren- 
der the BCGs 0.5 - 1 mag bluer than the g-r red-sequence 
(Bower et al., 1992, Baldry et al, 2004). This may impact 
the creation of large optical cluster catalogue based on the 
presence of a well defined red-sequence as well as on the 
presence of a BCG. 

Szabo et al. (2010) have created the largest available 
catalogue of optically selected clusters from the SDSS DR6. 
This cluster catalogue is based on the matched filter method 
(Dong et al., 2008) and does not include any colour selection 
of member galaxies. An interesting by-product of the clus- 
ter finder algorithm is the creation of the largest available 
catalogue of Brightest Cluster Galaxies with homogeneous 
photometry (and photometric redshift) without any selec- 
tion on colours. In the following we will refer to the Szabo 
et al. (2010) cluster catalogue and the related BCG cata- 
logue as synonyms, since each BCG is uniquely associated 
to a cluster. 

The main aim of this paper is to characterize such a 
catalogue by describing some tests done to check the its ac- 
curacy and characterizing the BCG properties in the colour- 
colour and colour-magnitude spaces. By means of a redshift- 
independent definition of blue BCGs, we will be able to 
quantify the bias that affects catalogues created via a colour- 
based selection. 

Moreover, with such a BCG catalogue we are in the 
position of pursuing several main goals. For instance, in 
this paper we extend the catalogue by including UV-optical 
colours of the BCGs thanks to Galex public data. By means 
of a positional cross-matching of our catalogue with pub- 
lished compilations of X-ray selected clusters we are in the 
position of extending Bildfell et al.(2008)'s and Pipino et 
al. (2009a) 's results to a larger sample of clusters/BCGs. 

A quantification of how cooling flows and blue BCGs 
impact cluster detection from Sunyaev-Zeldovich surveys is 
addressed in a companion paper (Pipino & Pierpaoli, 2010). 

The scheme of the paper is thus the following: in Sec. [2] 
we briefly describe the Szabo et al. cluster catalogue, briefly 
summarizing the BCGs selection and properties. In Sees. [3] 
4 and 5 we will present the characteristics of the BCG sam- 
ple in terms of luminosity functions, colours and redshift 
evolution of these properties as well as we compare them 
to other existing catalogues. We apply the catalogue to the 
study of the UV-optical colours and the X-ray properties of 
our BCGs in Sec. [6] Conclusions will be drawn in Sec. [7] 



2 THE CATALOGUE 

In this section we briefly summarize how the Szabo et al. 
cluster catalog that we use has been built. A detailed de- 
scription can be found elsewhere (Szabo et al., 2010). 
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Figure 1. Rodshift distribution of all the BCGs and those in rich 
clusters. 

2.1 Data 

The data on which Szabo et al. catalogue has been built 
are from SDSS DR6 (Adelman-McCarthy et al., 2008). All 
galaxy measurements were extracted from the Galaxy view 
on the CasJobs DR6 databas^]. Szabo et al. selected only 
galaxies that are detected in lxl binned images, that have a 
measurable profile, that are not saturated or contain peaks 
other than that provided by the estimator of the SDSS 
pipeline. 

The adopted photometric redshifts are based on a neu- 
ral network cc2 estimator and are available in the table Pho- 
toz2 (Oyaizu et al., 2008). We made this choice because of 
their more reliable error estimates and lack of evident colour 
biases as opposed to those available from Photoz (Csabai et 
al., 2008). For more details we refer to Szabo et al. (2010) 's 
paper. Absolute magnitudes are calculated by means of the 
kcorrect code v4.1.4 (Blanton & Roweis, 2007). In the fol- 
lowing we will always use SDSS model magnitudes. 

2.2 The cluster catalogue 

The matched filter method (Kepner et al., 1999) that Szabo 
et al. (2010) use is presented in detail by Dong et al (2008). 
In practice, it is a likelihood method which identifies clusters 
by convolving the optical galaxy survey with a set of filters 
based on a modeling of the cluster and field galaxy distri- 
butions. A cluster radial surface density profile, a galaxy lu- 
minosity function, and redshift information (when available) 
are used to construct filters in position, magnitude, and red- 
shift space, from which a cluster likelihood map is generated. 
The peaks in the map thus correspond to candidate cluster 
centers where the matches between the survey data and the 
cluster filters are optimized. The algorithm automatically 
provides the probability for the detection, best-fit estimates 

1 http://cas.sdss. org - /dr6 



of cluster properties including redshift, radius and richness, 
as well as membership assessment for each galaxy. Usage 
of the apparent magnitudes and the redshift estimates in- 
stead of simply searching for projected galaxy over-densities 
suppresses the foreground-background contamination. 

The cluster catalog is constructed with an iterative pro- 
cedure. The process starts from a density model of a smooth 
background with no clusters. For each galaxy position, we 
then evaluate the likelihood increment we would obtain by 
assuming that there is in fact a cluster centered on that 
galaxy. At each iteration, the cluster candidate which re- 
sulted in the greatest likelihood increase is retained. A list 
of cluster candidates then becomes available in decreasing 
order of detection likelihoods. The cluster richness A200 is 
then defined to be the total luminosity in units of L* inside 
r2oo, namely the radius inside which the mass over-density 
is 200 times the critical density. 

For the magnitude filter, Szabo et al adopt a luminos- 
ity profile described by a central galaxy plus a standard 
Schechter luminosity function (Schechter, 1976). 

Dong et al. (2008) showed that the selected cluster sam- 
ple is ~ 85% complete and over 90% pure for systems more 
massive than 1.0 x 10 14 /i _1 M (A200 ~ 50) with redshifts 
in the range 0.1-0.4. In order to have a reliable assessment 
of the BCG colours, we restrict the sample to galaxies in 
the redshift range [0.1,0.3], where the intrinsic scatter in the 
observed g-r colour is the smallest. The estimated cluster 
redshifts derived from maximum likelihood analysis show 
small errors with Az < 0.01. 

The final Szabo et al.'s catalogues has more than 69000 
entries. We refer the reader to the main catalogue paper for 
details on the cluster characteristics and comparison with 
other automated catalogues. A complete version of the cata- 
logue with the three brightest galaxies from which we derive 
the sample discussed in this paper can be found Szabo et al. 
(2010). 



2.3 A test of the accuracy of the catalogues 

We tested the accuracy of our selection by cross-matching 
the position of such an extended BCG sample with known 
coordinates of well studied BCGs taken from the literature 
(Crawford et al., 1999, Bildfell et al., 2008, Loubser et al., 
2008 as well as the SIMBAD database) whose position is 
within the region of sky covered by SDSS DR6 and that lie 
at 0.1 < z < 0.3. While we will mostly focus on BCG as the 
"brightest" member in every group or cluster of galaxies, it 
is important to consider also the 2nd and 3rd most luminous 
galaxies in a cluster (by using the r-band magnitude lumi- 
nosities) in the following exercise. This is needed because 
the above mentioned works either defined the BCG as the 
brightest galaxy in a band other than the r one or by their 
position in the cluster. Therefore there are a few cases in 
which their BCG is either the second or the third brightest 
in our definition. 

Our sample matches 73% of the BCGs in the Crawford 
et al. (1999) list with an accuracy of less than 0.3" in angular 
position and less of 0.01 in redshift space, and more than 
80% with if we require an accuracy of a few arc seconds, 
which is the typical error in the coordinates. Also differences 
in coordinates among authors (and the SIMBAD database) 
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for the same galaxy for the same galaxy may amount to a 
few arc seconds. 

In the remaining 20% of the cases: 

i) the known BCG that we want to match lies too close to 
the edge of the sky region covered by the SDSS DR6 and 
the AMF finder has problems in identifying a cluster there. 

ii) it is very close to the limit redshift range considered in 
this work. As an example we mention the case when the 
BCG to be matched has redshift 0.1, but its host cluster 
has a redshift < 0.1 and falls in the region where our cata- 
logue is not complete (and thus missed). 

iii) it is part of a sub-structure of a bigger cluster that it is 
not resolved by our cluster finder. 

A further complication is that, in the above mentioned lit- 
erature compilations, the authors often decide which galaxy 
to label as the BCG when galaxies with very similar lumi- 
nosity were present. In some cases their selection has been 
done on the basis of the presence of an extended stellar 
halo which clearly made the BCG candidate a cD galaxy. 
In other cases, the galaxy closer to the X-ray centre of the 
cluster has been selected. Finally, there are cases in which a 
galaxy cluster is made of sub-clusters in the act of merging, 
each one with its own BCG. Therefore it is not surprising 
that we have a fraction of cases in which the positional best 
match with a known BCG has been obtained by the second 
or the third most luminous galaxies in the entire Szabo et 
al.'s BCG sample. Also, we have cases in which the first and 
the second brightest galaxies of a given cluster in our entire 
BCG sample match the BCGs of two known (merging) sub- 
clusters. A more detailed assessment of these cases will be 
the topic of a forthcoming study. 

To our knowledge, the Szabo et al. cluster catalogue is 
the only one tested in comparison to known cluster and BCG 
positions. Other works (e.g. Koester et al., 2007) limited 
their analysis to matching the position of their clusters with 
the centre of a known cluster. Here we show that our catalog 
includes the single galaxies classified as BCGs from previous 
works. 

2.4 The final BCG sample 

We will focus on BCG as the brightest member in every 
group or cluster of galaxies in the r band. We further exclude 
from the analysis galaxie s whose error on the g-r colour, 
expressed as 3 • \/og~+~(j?, exceeds 0.3 mag. They represent 
about less than 1% of the galaxies in the redshift range 0.1- 
0.3. The final sample comprises of 14344 galaxies. 

The number of BCGs (Fig. [TJ increases with redshift 
as the underlying cluster redshift distribution. The richness 
of the cluster does not play a role: the histogram for rich 
clusters tracks the one for the whole population. 



3 BCG CHARACTERIZATION 

In this section we will focus on the BCG characterization in 
terms of luminosity, colours and explain how the inclusion of 
blue BCGs make Szabo et al. (2010) catalogue different from 

2 as opposed to other selection criteria such as being the highest 
likelihood member or the closest to the estimated cluster centre. 
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Figure 2. Distribution of galaxies in the plane r-band magnitude 
versus photometric redshift: all BCG (points and solid contours, 
upper panel) and galaxies in clusters with richness larger than 
50 (solid contours, lower panel). BCGs bluer than 0.3 mag from 
the red-sequence at their redshift arc shown by dotted iso-deusity 
(number of galaxy/bin area) contours. Dashed lines: lcr region 
around the mean relation by Loh & Strauss (2006) for BCGs in 
SDSS. Compare with Koester et al (2007) Fig.3. 

others in the literature. For a comparison of BCG luminosity 
distributions between catalogues we refer to Szabo et al. 
(2010). From Sec. 13. 21 we anticipate that we refer to galaxies 
that are 0.3 mag below (i.e. bluer) than the g-r red-sequence 
at their respective redshifts as blue BCGs. 



3.1 BCG luminosity 

In Fig. [2l we compare our BCGs to the BCG magnitude- 
redshift relation inferred by Loh & Strauss (2006) for LRG 
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in the redshift range 0.12-0.38. The dashed lines bracket the 
la scatter around this relation. At high richnesses (lower 
panel) the spread in our BCG sample is comparable with a 
~3tr scatter around Loh & Strauss' relation as well as the 
mean trends look very similar to each other. In particular, 
we find that r ~ 12.6 ■ z + 14.7. A substantial population 
of blue BCGs is present when considering the entire sample 
(upper panel) and offsets our distribution towards fainter 
magnitudes with respect to Loh & Strauss (2006) findings. 

In Fig. [3] (upper panel) we show the distribution in lu- 
minosity in the r-band for the BCGs in the redshift range 
[0.1,0.3] as a thick solid line. This curve slightly deviates 
from a Gaussian distribution with mean equal to the aver- 
age M r in the same redshift range and a ~0.5 mag which 
is commonly adopted as the BCG luminosity function (e.g. 
Hansen et al., 2005). As we will discuss in more details later 
in the paper, bluer BCG tend to populate poorer systems. 
Since the BCG luminosity scales with the cluster richness 
(Lin & Mohr, 2004, Hansen et al., 2009), blue BCGs (dashed 
line) tend to be fainter than the average. Note that here we 
are showing the distribution for the entire BCG sample, ir- 
respective of the galaxy redshifts. In fact, when we look at 
the distributions in two redshift bins for the whole sample 
(dash-dotted lines), we notice that the distributions are nar- 
rower. 

The average M r increases (conversely the luminosity 
decreases) at smaller redshift. In particular we find that 
the mean M r scales linearly with log (time) as expected 
from pure passive evolution (e.g. Tinsley, 1980, Nelson et 
al., 2001). The average luminosity drops by ~20% from red- 
shift 0.3 to 0.1. Also highlighted (lower panel) are the distri- 
bution functions in clusters with richness above 50. A com- 
parison between the two panels makes evident that the tail 
at faint magnitudes is due to BCGs in low richness clus- 
ters. While a study of the BCG luminosity function clearly 
deserves more attention, we can conclude that we found a 
distribution that it is broadly consistent with a Gaussian 
shape and featuring a redshift evolution in agreement with 
a passive evolution of the BCG population (e.g. Brough et 
al., 2002). A more quantitative assessment of the BCG lu- 
minosity function and its evolution would require a more 
careful treatment of the data. Here we rely on the absolute 
magnitude in the r-band calculated by means of the tem- 
plate fitting approach (Blanton & Roweis, 2007, Csabai et 
al. 2008) but with photometric redshift derived from a neu- 
ral network estimator (Oyaizu et al., 2008), namely not is 
self-consistent way. Furthermore no evolutionary corrections 
have been applied. Moreover, it is known that sky subtrac- 
tion errors in the SDSS pipeline may significantly affect the 
magnitude of the brightest objects (Adelman-McCarthy et 
al. 2008, and references therein). Finally the catalogue is not 
complete at the poor richness end, and this implies that the 
low-luminosity tail of the luminosity distribution is not cor- 
rectly represented. Finally, the effect of small number statis- 
tics is clear in the distribution for rich (A200 > 50) cluster. 

In concluding the section, we add that Fig. [3] suggests 
a cut at M r < —22.5 if one wants to use a sample almost 
made by red BCGs. 
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Figure 3. Upper panel. Solid thick line: luminosity distribution 
the r-band for BCGs in the redshift range [0.1,0.3]. Lighter thick 
line: Gaussian distribution with mean equal to the average of the 
above distribution and <r=0.5 mag. The function for blue BCGs 
(dashed line) and the subsets of all BCGs divided into two redshift 
bins (dot-dashed lines) are also shown. Lower panel: As above, 
but only for rich clusters (A200 > 50). 



3.2 Colour classification 

In order to understand the difference between a cluster (and 
BCG) catalogue which relies on the matched filter method 
with respect to a colour classification, it is important to 
define and characterize the fraction of blue BCGs. In this 
work, galaxies 0.3 mag below (i.e. bluer) than the average 
g-r colour (observed frame) for early types at their respective 
redshifts are considered as blue BCGs. Since these galaxies 
obey to the color-magnitude sequence, and since we focus on 
a narrow range in luminosities (at the high mass/luminosity- 
end), in practice our requirement implies that these blue 
BCGs will lie below the red-sequence at their given redshift. 
Indeed, we chose such a cutoff since 0.3 mag is more than 6cr 
off the mean value of the color-magnitude relation (Bower 
et al., 1992). In practice, we self-consistently evaluate the 
average g-r at a given redshift by means of our BCG sam- 
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(brighter points) than our cutoff, which is roughiy given by 
the reiation g — r ~ 2.86 ■ z + 0.35 (solid line) is larger than 
that of the galaxies on the red side of the average sequence 
(g — r ~ 2.86 ■ z + 0.65, dashed lines) . 

In order to guide the eye we also plot - bracketed by 
lighter dashed lines - the la region around the mean relation 
found by Koester et al. (2007) for their BCGs (c.f.their Fig 2) 
which follow more closely the mean red-sequence evolution 
as a function of redshift. On the other hand, the bulk of 
our BCGs in clusters of richness above 50 matches very well 
the relation for red and dead early type galaxies. A number 
of blue outliers are still present though. We further discuss 
these galaxies below. 
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Figure 4. g-r color as a function of photometric redshift for both 
all BCGs and those in rich clusters. BCGs bluer than 0.3 mag 
from the red sequence at their redshift (i.e. below the solid line, 
in the "blue" region) are presented by means of brighter points. 
We also highlight them using the asterisks symbol in the panel for 
high richness clusters, whereas squares highlight spirals. Dashed 
lines: la region around the mean relation found by Koester et al. 
(2007) for their BCGs (their Fig 2, see also Blanton & Roweis, 
2007). 



pie and we use it as a zero point for measuring the offset 
(i.e. the blueness) of the single galaxies. Not surprisingly, it 
turns out that this coincides with the mean trend of the g-r 
versus redshift curve reported by Blanton & Roweis (2007, 
c.f. their Fig. 3), hence basically with the locus where galax- 
ies in the Cut I sub-sample of the Luminous Red Galaxy 
(LRG) sample (Eisenstein et al., 2001) lie. In other words, 
it is very likely that LRGs that are also BCGs in their own 
cluster, end up as BCGs in the Szabo et al. catalogue. The 
converse is not true, since the colour cuts in the LRG sample 
get rid of the galaxies that, albeit very bright, are 0.3 mag 
bluer than the average g-r colour. The fraction of the galax- 
ies bluer than the average with respect to the total number 
of BCGs gives us an estimate of the colour bias - possibly 
introduced by recent star formation (and cooling flows) - 
into red-sequence based cluster catalogues. 

In Fig. [4] we show the g-r color as a function of pho- 
tometric redshift for all BGCs (upper panel), whereas we 
focus on those in rich cluster in the lower panel. We note 
a clear asymmetry in that the number of BCGs 0.3 bluer 



3.2.1 The blue fraction 

We now quantify the fraction of blue galaxies in terms of 
redshift, richness and dominance. In Fig. [S] we show the dis- 
tribution in the offset in magnitudes from the average g-r in 
different redshift slices. Whilst the large majority of BCGs 
display a negligible offset and follow the same trend of the 
luminous red galaxies studied by Blanton & Roweis (2007), 
the histograms show a clear tail on the blue side (positive 
offset in our figures). As expected from Fig. 2] there is a 
clear asymmetry in the number of BCGs bluer than the 
red-sequence, with respect to those redder that we showed 
in Fig. 2] Note, however, that the curves become more sym- 
metric and the relevance of the blue BCGs diminishes if we 
take into account only rich clusters and as we move to lower 
redshift. In particular, we find that the overall fraction of 1st 
ranked galaxies bluer (redder) than 0.3 mag is 8.6% (0.2%). 
Such a fraction does not change if we restrict ourselves to 
clusters with richness above 30, and decreases to the 6% 
above a richness of 50. In terms of redshift bins, the blue 
fraction goes from 5% in the redshift range 0.1-0.2 to 10% 
in the redshift bin 0.2-0.3. Thus blue BCGs tend to popu- 
late poor clusters and their fraction slightly increases with 
redshift. Note that at higher redshifts both the scatter in 
the g-r colour and the its error tend to increase. 

The overall blue fraction is lower than the 25% of BCGs 
with blue optical cores found by Bildfell et al. (2008), be- 
cause in not all cases the spatial extent of the blue core is 
large enough to make the entire galaxy blue0. The fraction 
of 1st ranked galaxies bluer than 0.5 mag is 2.8%. 

A common colour-cut used by discriminate between 
early type morphologies (which BCGs belong to) from later 
type at u-r>2.22 has been derived by Strateva et al. (2001) 
on an earlier release of the SDSS. This requirement is satis- 
fied by the 95% of our BCGs in clusters with richness above 
30, in agreement with Strateva et al. (2001), who found that 
97.6% of their (spectroscopically classified) early type galax- 
ies were above the u-r=2.2 threshold. On the other hand, 
only the 13% of those that are blue according to our clas- 
sification make this threshold. In agreement with us, Choi 
et al. (2007), found that about 10% of early type galaxies 
with evidences of a blue optical core in a volume limited 
sample of SDSS has a u-r colour bluer than 2.22. Most of 
these galaxies live in low-density environment as we found. 



3 There are remarkable cases (e.g. Abell 1835) where the BCG 
can be 0.5-1 mag below the red-sequence 
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Figure 5. Offset (mag) from g-r color-magnitude relation as a 
function of redshift and offset distribution in redshift slices. The 
shaded area emphasizes the tail of blue BCGs. 

Unfortunately a cross-match with their catalogue is impos- 
sible because it does not overlap in redshift with ours, being 
theirs at redshift below 0.1. 

The presence of g-r blue BCGs can be inferred by us- 
ing other colours. For instance, galaxies blue in g-r have also 
quite blue u-g colours, similar to those of star forming galax- 
ies. This is shown in Fig. [6] In this plot the solid diagonal 
line divides the plane at u — r — 2.2: the upper right portion 
is the typical locus of LRG. The solid contours, that pertain 
to all the BCGs, peak in this region, presenting, however, 
a large scatter. One natural reason for this is the presence 
of the blue BCGs (note secondary peaks at g-r~u-r~l that 
coincide with the peak of the blue BCG distribution (dot- 
ted contours). Other reasons might be traced back to the 
fact that there are observational errors (not taken into ac- 
count in this simple analysis). Moreover galaxies are made 
by mixtures of stars of different ages and metallicities. This 
creates a scatter in the colours of otherwise similar galax- 
ies (for instance galaxies with the same mass and about the 
same age, see e.g. the Maraston et al., 2009, but for the LRG 
case). Moreover, galaxies are still likely to have a spread in 
mean ages (namely, it is very unlikely that they have formed 
simultaneously at, say, z=3). Blue BCGs cluster in a com- 
pletely different region of this diagram. A colour-cut at u- 
r=2.2 would get rid of most of such galaxies, thus biasing 
the study of the BCG population as a whole. Other colours, 
such as the i-r, are less useful, since the regions occupied by 
blue and red galaxies tend to significantly overlap. 

3.2.2 The blue fraction in galaxies with spectra 

The fraction of blue BCGs in the whole sample goes down 
to 3.7% when considering only the galaxies that have a spec- 
troscopic redshift. In particular, the fraction of blue BCGs 
in cluster with richness above 50 that have a spectroscopic 
redshift decreases to the 1%. In this latter case, all the galax- 



Figure 6. Distribution of galaxies in the (g-r) - (u-g) plane. The 
solid diagonal line divides the plane at u — r = 2.22: the upper 
right portion is the typical locus of LRG. The solid contours per- 
tain to the entire sample of BCGs. The blue BCGs are presented 
by the small points and highlighted by dotted contours. 

ies with g — r < 0.7 would disappear from the lower panel in 
Fig. [Hand their fraction would be strongly suppressed in the 
upper panel. We cannot exclude that spectra were preferen- 
tially taken for galaxies sitting on the red-sequence, hence 
artificially removing blue galaxies, since we have spectro- 
scopic redshifts for just one third (5642/14344) of the BCGs 
in the range studied. However, a close inspection to galax- 
ies in rich clusters tells us that, while the difference between 
spectroscopic and photometric redshift tend to be very small 
(< 0.02) in red galaxies, it is not uncommon to find "blue" 
galaxies whose spectroscopic redshift is 0.1-0.2 lower than 
the photometric estimate. Given the quite steep variation of 
(g-r) with redshift (c.f. Fig. [4]), such a galaxy might be per- 
fectly on the red-sequence when the correct (spectroscopic) 
redshift is used. Therefore, even considering the overall ac- 
curacy of state-of-the-art photometric redshifts estimates, 
we cannot exclude that a fraction of the blue BCGs are so 
because of the specific error in their photometric redshift. 

Therefore, a conservative final estimate of the blue frac- 
tion of BCGs in the redshift range 0.1-0.3, at richnesses 
above 50, is between 1 and 6%. Such estimate is in agree- 
ment with the results by Edwards et al. (2007), who selected 
BCGs by means of K-band magnitudes and found that ~3% 
of those had bluer (u-r) colour than expected. 

3.2.3 The blue fraction at higher redshift 

Despite the focus on the present paper is on galaxies in the 
redshift range 0.1 — 0.3, the parent BCG catalogue has a 
wealth of data for higher redshift objects. In this section 
we briefly quantify the blue fraction in such a regime. The 
reader should note that the colour selection at z ~ 0.3 is 
somehow more complicated. In the first place, the complete- 
ness of the catalogue quickly decreases at z > 0.4 , therefore 
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we cannot properly assess the fraction of blue BCGs. Fur- 
thermore, the the g-r - redshift curve flattens out already 
at z ~ 0.35. Both intrinsic spread and errors in the colours, 
as well as the possible presence of evolutionary effects on 
the color-magnitude relation, enhance the fraction of blue 
BCGs at higher redshift. However, taken at face value, our 
results in the redshift range 0.3 — 0.43 would indicate that 
the overall blue fraction increases to 18% (14% for clusters 
with A200 > 50). A more conservative estimate can be ob- 
tained if we consider the fraction of galaxies ~ 5cr below the 
mean g-r at z~0.4 as the blue ones, namely 0.5 mag instead 
of 0.3. In such a case the blue fraction is 11% (8.5% clus- 
ters with A200 > 50). These values are very close to those 
attained in the 0.1-0.3 redshift range. A similar increase in 
the fractions of BCGs with optically blue cores and emission 
lines has been observed by Wang et al. (2010) and Crawford 
et al. (1999), respectively. As we will see in the following 
section, our results confirm and extend the correlation be- 
tween BCG colour and cluster mass (e.g. More et al., 2010, 
Loh et al., 2010) at redshifts z>0.1. 

3.3 Morphologies 

Until this point we did not take into account the morphology 
of the BCGs. A more detailed study on the morphologies in 
the Szabo et al. catalogue is underway (MacKenzie et al., in 
prep). It will tell us whether the blue BCGs are truly early 
type galaxies or if the sample is contaminated by late type 
(i.e. bluer star forming) galaxies. These may be identified as 
BCGs if sufficiently bright members of clusters lack, e.g., a 
well defined red sequence or a prominent early type central 
galaxy. In particular, for application to galaxy formation 
studies, it can be worthwhile investigating if the blue frac- 
tion can be associated to other galaxy types as opposed to 
elliptical ones. For instance, lenticular galaxies are BCGs in 
known clusters, while spiral galaxies may dominate groupfl 
Finally the blue colours may reflect the star formation in 
interacting systems. Here we mention some preliminary re- 
sults: in the richest (A200 > 100) clusters of our sample, 
93% of the galaxies are single or interacting ellipticals, the 
rest being "uncertain" according to a visual morphological 
classification. In particular, the interacting ellipticals alone 
amount to 34%, similar to what found for z ~ massive el- 
lipticals by Kannappan et al. (2009). According to Kannap- 
pan et al. (2009), these galaxies populate low-to-moderate 
density environments. The distributions of these interact- 
ing early type galaxies around the average g — r colour at 
a given redshift and in absolute magnitude is very similar 
to the one that single ellipticals exhibit. For the above rea- 
sons, we suggest a cut at A200 > 100 for studies focussing 
on red early-type BCGs in massive galaxy clusters. A vi- 
sual/spectral classification of the blue BCGs in A200 > 50 
clusters of our sample, tells us that the percentage of ellipti- 
cals decreases to 70%. 5% of these blue galaxies are indeed 

4 By adding galaxies in this redshift range, the total number of 
BCGs would nearly double. 

5 It is worth reminding the reader that relation between A200 
and the cluster mass has a substantial scatter below A200 = 50 
(Fig. 6 in Dong et al., 2008), such that ~ 2 ■ lO 14 h -1 Af0 poor 
clusters and ~ 5 ■ 10 ls h~ Mq groups can be assigned the same 
richness. 



spirals typically associated to A200 ~ 50 clusters (squares 
in Fig. [4|. Being associated to poorer systems, these spi- 
rals are on average fainter. Therefore we suggest a cut at 
M r < —22.5 if one wants to use a sample almost made by 
BCGs that are early-type galaxies. The morphology-richness 
trend seems to follow the increase in the blue fraction at 
lower richness that we discussed above. It is worth noting 
that similar findings are also reported by other studies of the 
correlation between BCG properties and cluster mass. For 
instance, More et al. (2010), using the kinematics of satellite 
galaxies to infer the halo mass, found that, at a given galaxy 
luminosity, red central galaxies tend to occupy more mas- 
sive haloes than the blue ones. Similar results are reported 
by Loh et al. (2010) by means of the two-point correlation 
function as well as the NUV-r colours (see below). Hence, 
they are able to verify that not only red galaxies but also 
objects in transition between the blue and the red sequence, 
preferentially reside in more massive haloes. These indepen- 
dent results have been obtained below z ~ 0.1, therefore we 
can confirm the reported trends and extend their validity 
out to z = 0.3. 

Finally, we cannot exclude that some objects in the en- 
tire sample are background/foreground field object with pe- 
culiar colours. 



4 THE BCGS AS THE THREE BRIGHTEST 
GALAXIES IN A CLUSTER 

In this section we consider the three richest galaxies as the 
potential BCGs of a given clusters, because in know local 
clusters the brightest is not always the dominant one. More- 
over, massive clusters often show signs of on-going merg- 
ers and clear substructures, each one with its own BCGs. 
Finally the importance of comparing the properties of the 
brightest galaxy to the second and the third brightest ones 
is related to potential applications of the catalogue in order 
to constrain galaxy formation models. 



4.1 The role of the dominance 

We first turn our attention to the role of the dominance of 
the BCG. We express it in terms of the difference in r-band 
apparent magnitude V2nd — ri s t between the 2nd and the 
1st ranked galaxies. The average dominance in the redshift 
range 0.1-0.2 is 0.5 mag in the r-band, a value within the 
range given by Loh & Strauss (2006) 0. Also our lo" disper- 
sion and the mild decrease of the dominance with redshift 
are in agreement with their findings. In Fig. we show the 
distribution of galaxies in the cluster richness - dominance 
plane. As expected from other studies (e.g. Loh & Strauss, 
2006), poorer systems exhibit a larger average dominance 
than richer clusters. In the former systems the BCG, albeit 
less luminous, contributes a larger fraction of the total light 
than in the latter. In rich clusters, instead, the fraction of 
the total light coming from the BCG is much smaller and 
such a difference can tell us something on the path that lead 



6 Note that a direct comparison cannot be made, since they fo- 
cused only on the LRGs 
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Figure 7. Distribution of galaxies in the cluster richness versus 
1st ranked BCG dominance plane. 
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Figure 8. Distribution of galaxies in the offset from the red- 
sequence versus BCG dominance plane. The blue BCGs are pre- 
sented with dotted contours. The numbers give the actual number 
of galaxies within that contour. 

to the formation of the BCG (e.g. Lin & Mohr, 2004). In par- 
ticular, the dominance as a function of richness, as well as 
the luminosity of the brightest and 2nd brightest galaxies as 
a function of the dominance are very promising tools to con- 
strain galaxy formation models and to discriminate among 
apparently similar recipes for the BCG star formation his- 
tory (Smith et al., 2010). 

Finally, in Fig. [8] we investigate the effect of the domi- 
nance on the 1st ranked galaxy blueness, quantified in terms 
of the offset in magnitude from the red sequence at the 



galaxy redshift. No significant trends are found, except for 
the fact that cluster with low dominance are much more 
abundant than clusters where rmd — r\ a t exceeds 1 . 

4.2 The blue fraction 

As far as the blue fractions under such a broader BCG defi- 
nition are concerned, we find that it is 12% overall, 9.8% in 
clusters with A200 > 30 and 8.2% in clusters with A200 > 50. 
Such values are ~1.5 times larger than the fractions at- 
tained when considering only the brightest cluster galaxy. 
This fact is a consequence of the colour-magnitude relation, 
namely, the fact that in a given cluster the brightest early 
type galaxy is also the reddest. Therefore, in a given cluster, 
lower luminosity galaxies are likely to be bluer even if they 
have an early-type morphology and no signs of star forma- 
tion. However we cannot exclude a slightly higher presence 
of contaminants as bright spiral galaxies that are in the pro- 
cess of being accreted in the cluster. The trend of the blue 
fraction with redshift in clusters with A200 > 50 for such 
a broader definition of BCGs mirrors the one found when 
considering only one (the brightest) galaxy per cluster. 



5 A COMPARISON WITH MAXBCG BCGS 

The presence of blue BCGs is an important feature of the Sz- 
abo et al. cluster catalogue. Hence, it is important to present 
a quantitative comparison with a widely used catalogue as 
maxBCG (Koester et al., 2007) based on a strict colour se- 
lection. We briefly remind here that our BCG is defined as 
the brightest galaxy in the r-band that likely belongs to a 
cluster, even though it does not necessarily sit at the cluster 
centre. On the other hand the maxBCG algorithm requires 
a recQ and bright galaxy and at least other ten red and less 
luminous galaxies within ~ 1 Mpc to identify a cluster. The 
seed galaxy is hence the BCG and the central galaxy at the 
same time. The cross-matching between the two catalogues 
has been done by searching for maxBCG BCGs that are 
also one of three brightest galaxies in one of our clusters. 
In particular, we find that more than 4300 maxBCG are in 
the entire Szabo et al.'s BCG sample. In the majority of the 
cases, there is a one-to-one correspondence, in the sense that 
only one maxBCG BCG belongs to a Szabo et al.'s cluster. 
We define these clusters as those that match the maxBCG 
catalogue. For a more thorough description of the matching 
procedure and a complete analysis of the outcome we refer 
to Szabo et al.'s paper. 

In the upper panel of Fig. [9] we show the distribution 
in the offset in magnitudes from the average g-r colour for 
all our BCGs (solid line) and those that are in clusters in 
common with maxBCQj (dashed line) as well as those that 
are not in common (dotted line). Note that the comparison 
has been made by taking only Szabo et al's BCG in the sky 

7 We refer to Koester et al., 2007 for the actual colour cuts; here 
it suffices to say that the criterion is such that the galaxy colours 
are within the dashed lines of Fig. 4. 

8 Note that the fact that a fraction of our clusters matches 
maxBCG clusters does not imply that our 1st ranked BCG in 
these clusters always coincides with the BCG of the maxBCG 
catalogue. 
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Figure 9. Upper panel: Offset (mag) from g-r color-magnitude re- 
lation for all our BCGs (solid line) and those that are in clusters in 
common with maxBCG (dashed - see text) and those that are not 
(dotted line). Lower panel: Offset (mag) from g-r color-magnitude 
relation for all our BCGs (solid line, as in the other panel) versus 
the offset from the color-magnitude relation of maxBCG BCGs 
that belong to our clusters but that do not coincide with our BCG 
(dashed line). 



regions where SDSS DR6 and DR5 (on which maxBCG is 
based) overlap and in the redshift range 0.1-0.3. We dis- 
play normalized histograms in order to emphasize the tail 
at positive (blue) offsets in our BCGs. As expected from 
the maxBCG algorithm, the colours of maxBCG BCGs that 
also are BCGs for our clusters do not differ from those of 
maxBCG BCGs that belong to our clusters without being 
the brightest member. The distribution of our BCGs, in- 
stead, shows a more evident tail, because clusters that do 
not have a maxBCG counterpart tend to be poor systems 
(see Szabo et al., 2010), where the fraction of blue BCG is 
somehow larger (see above). The offset distribution for the 
galaxy that Szabo et al. classify as 1st ranked BCG in a 
given cluster and that for the BCGs selected according to 
the maxBCG method, basically coincide but for the tail with 
blue galaxies. Therefore, in clusters that Szabo et al and 
maxBCG share, the BCG selection is almost identical. In 
7% of the cases maxBCG miss the brightest object because 
is more than 0.3 mag away from the red-sequence. However, 
differences are found already at 0.2 mag away from the red 
sequence. Such a fraction (7%) is close to 6%, namely is 
the fraction of blue BCG in clusters with richness > 50 (see 
above) , and it is a sensible value because most of our matches 
with maxBCG clusters are in the high richness regime (see 
Szabo et al., 2010). Since we know that the galaxy classified 
as BCG by Koester et al. (2007) is still one of the brightest 



(but not our 1st ranked) galaxies in our clusters, this means 
that maxBCG uses the 2nd or the 3rd brightest galaxy as 
seed for their clusters. Indeed, in the lower panel of Fig.|5]we 
plot the offset distribution for all our BCGs (solid line, as in 
the other panel) versus the offset from the color-magnitude 
relation of maxBCG BCGs that belong to our clusters but 
that do not coincide with our BCG (dashed line). 

In conclusion, if more than 80% of maxBCG rich clus- 
ters have a counterpart in our cluster catalogue (see Szabo 
et al., 2010), we can confirm that they miss a fraction of 
BCGs bluer than expected (namely 0.1 - 0.2 mag off the red- 
sequence, see Koester et al., 2007). In other words, the colour 
criteria used to select maxBCG BCGs may not bias cluster 
catalogues but certainly affect BCG catalogues to the 10% 
level. A caveat is that, in order to make the comparison on a 
the same "frame" , we used colours and " model" magnitudes 
for maxBCG BCG as provided by the SDSS DR6. Since the 
original maxBCG catalogue has been built from the SDSS 
DR5, used "cmodel" magnitudes and different photomet- 
ric redshift estimator (and hence different k-corrections) the 
exercise presented here does not provide a colour character- 
ization of maxBCG BCGs in a strict sense. 

On the basis of the maxBCG BCG colour alone, instead, 
we cannot explain why a fair number of maxBCG clusters 
is not matched by us. As shown by Szabo et al., the dif- 
fering completeness above a given richness and the differing 
definition richness between the two catalogues hamper a 1:1 
matching. Also, we cannot exclude that we include groups 
where the BCGs are spirals. 



6 APPLICATIONS OF THE CATALOGUE 

After presenting the general characteristics in terms of op- 
tical colours of the Szabo et al. BCG catalogue, we discuss 
some possible applications. To this aim, we also augment the 
available data per each cluster/BCG by cross-matching with 
catalogues at other wavelengths. In particular, we briefly 
present some results on the UV-optical colours of our BCGs 
and the X-ray properties of their host clusters in connec- 
tion to the galactic blueness. In these cases, since we are 
matching our sample with existing databases that are nei- 
ther complete nor homogeneous, and whose sky coverage 
only partially overlaps with ours, we will not yield a com- 
plete catalogue. The following examples cannot, therefore, 
be taken as a characterization of our entire catalogue. Hence, 
in this paper we use them just to highlight the potential ap- 
plication of our BCG sample in spectral regimes other than 
the optical. 



6.1 BCGs in the UV 

In contrast to the optical spectral range, the UV is highly 
sensitive to even small fractions of young stars (younger than 
about a Gyr), making it an excellent probe of the low-level 
recent star formation (RSF) that is expected in elliptical 
galaxies (i.e. star formation within the last Gyr, contributing 
up to a few percent of the stellar mass of the galaxy, e.g. 
Kaviraj et al., 2007). As explained earlier, we expect this 
to happen also in BCGs, being them the most massive early 
type galaxies, since a fraction of them has blue cores (Bildfell 
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et al, 2008, Rafferty et al., 2008) and shows emission lines 
(Crawford et al., 1999, Edwards et al., 2007, 2010). 

Unfortunately the SDSS database does not provide us 
with colour gradients to understand if a BCG in our sample 
has a blue core due to recent star formation. Nor can the 
RSF always lead to blue galaxies in terms of their g-r colour. 
Therefore in this section we will make use of the results put 
forward by Pipino et al. (2009a, see also Wang et al., 2010). 
They found that every BCG which has a blue UV-optical 
colour also shows a blue-core in its optical colour profile. 
Conversely, BCGs that lack blue cores and show monotonic 
colour gradients consistent with a decrease in metallicity 
with radius typical of old elliptical galaxies are red in the 
UV. Pipino et al. (2009a) interpreted this as evidence that 
the UV enhancement in the blue BCGs is driven by recent 
star formation and not from old evolved stellar populations 
such as horizontal branch stars. The recent star formation 
in the blue BCGs typically has an age less than 0.5 Gyrs 
and contributes mass fractions of less than a percent. In a 
sense, the fraction of blue BCG derived in the previous sec- 
tions is a lower limit on the RSF in BCGs (or an estimate 
of the fraction of BCGs with intense enough RSF to have 
their optical colours offset from the red-sequence). The frac- 
tion derived in this section will be, instead, related to the 
percentage of BCGs experiencing RSF at a very low level. 

We cross- matched the entire BCG sample with publicly 
available UV photometry from the GR4 and GR5 data re- 
lease of the GALEX mission (Martin et al. 2005) and we find 
a counterpart for our entire sample of BCGs in roughly one 
third of the cases. In particular, we retrieved data by means 
of the cross-matched Galex GR4+GR5 - SDSS catalogue 
available in the Galex archive. The positional matching, per- 
formed within 5", returns nearly 5000 objects. The remain- 
der of the galaxies are not observed partly because the Galex 
sky coverage only partially overlaps with the SDSS footprint 
and because most of the overlap area has been only observed 
with All sky Imaging Survey (AIS - short exposure times) 
rather than with the deeper surveys (Martin et al., 2005). 
We note that, in general, galaxies harboring recent star for- 
mation are more likely to be detected than non recent star 
formation-galaxies (since their UV flux will be higher), es- 
pecially at higher redshift. Thus at any redshift the chances 
of detecting a UV blue BCG is higher than its UV red coun- 
terparts. 

Fig. [10] shows that optically blue BCGs systematically 
have a UV excess (i.e. NUV-r ~ 41)- I n particular, we find 
that while only the 14.7% of BCGs (that have a GALEX 
counterpart) in cluster with richness higher than 30 exhibit 
a NUV-r colour below 3, this fraction rises up to 71% when 
considering only the optically blue ones. Very similar frac- 
tions (12.5% and 72%, respectively) are obtained consider- 
ing BCGs in clusters whose richness exceeds 50. Note that 
in the NUV-r vs r color-magnitude diagram the dispersion 
is much larger than typical optical colours colour-magnitude 
relation; therefore the threshold at NUV-r=3 used here em- 
phasizes the galaxies with the bluest UV-optical colour. 

Choi et al. (2007, c.f. their Fig. 6) show that colour gra- 



9 It is worth reminding that NUV-r <5.3 is a conservative limit 
that discriminates between UV excess caused by recent star for- 
mation, and that caused by old stars (Yi et al., 2005). 
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Figure 10. NUV-r colour versus offset in magnitudes from the g-r 
red-sequence (observed frame). Points: 1st ranked BCGs. Squares: 
blue BCGs. Note: to improve the quality, only a limited number 
of galaxies in the region with g-r~0 and NUV-r~6 is plotted. 
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Figure 11. NUV-r colour versus u-r colour. Symbols as in Fig. 1101 



dients correlate with u-r colour for SDSS early type. They 
find that ~ 10% of morphologically classified ellipticals have 
positive colour gradients and u-r < 2.2. In particular, they 
find gradients that are positive (blue cores) at u-r < 2.2 
and becoming more and more negative (red and dead ellip- 
ticals) as the u-r colours gets redder. This is exactly what 
we find and show in our Fig. [TTJ if we assume that there is 
a one-to-one correspondence between the optical blue core 
and the UV excess (Pipino et al., 2009a). This correspon- 
dence is suggestive, however a more detailed work is needed 
to establish it on solid grounds. For instance, one should 
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perform a deeper analysis that includes a visual morpholog- 
ical classification, an inspection of all the matches to ensure 
that we do not include UV light coming from close compan- 
ions, foreground and background (including lensed) objects. 
This is left to a forthcoming study (MacKenzie et al., in 
prep.). 

6.2 BCGs in the X-ray 

Bildfell et al. (2008) found that the presence of optical blue 
cores in 25% of its BCG sample is directly linked to the 
X-ray excess of the host clusters. Moreover the position of 
these BCGs coincides with the peak in X-ray emission. Their 
interpretation is that the recent star formation in BCGs is 
associated with the balance between heating and cooling in 
the ICM in the sense that the clusters that are actively cool- 
ing are forming stars in their BCGs. The aim of this section 
is to confirm the results by Bildfell et al. (2008) and show 
how our BCG catalogue can be successfully exploited for 
studies of X-ray clusters and intergalactic medium proper- 
ties. For this purpose, we use the entire catalogue of (1st 
ranked) BCGs from Szabo et al., without any redshift limit. 

We make use of the publicly available data made pos- 
sible by the ACCEPT project (Cavagnolo et al. 2009). The 
catalog comprises 239 galaxy clusters with accurate tem- 
perature, density, entropy and pressure profiles reduced in 
a homogeneous way from public Chandra data. The catalog 
covers the temperature range 1-20 keV with redshifts rang- 
ing from 0.05 to 0.89. We refer to Cavagnolo et al. (2009, 
and references therein) for details on the data reduction and 
further catalog specifics. Here we note that their catalog 
is neither flux limited, nor volume limited. The matching 
procedure is similar to what done in Sec. 12.31 however a 
more detailed scrutiny is required. In the first place, in or- 
der to maximize the number of matches we do not limit the 
analysis to BCGs in the redshift range 0.1-0.3. We discard 
matches that have a difference in redshift larger than 0.03 
if the galaxy spectroscopic redshift is available, 0.1 other- 
wiscFl. In this latter case, however, we further check that 
the galaxies are associated to the clusters. In particular, we 
visually inspect the SDSS images in order to study the po- 
sition of the galaxy in relation to the literature position of 
the X-ray cluster. In practice, the match is rejected when a 
clear visual over-density of galaxies is found around the lit- 
erature X-ray cluster position but the alleged BCG is instead 
isolated and more distant. 

Such a procedure returns 38 matches of which 35 within 
200 kpc from the cluster centres (Fig. I12|l . This is a well- 
known properties of BCGs, i.e. they are typically located 
within few arc seconds from the cluster centres. 

We compute the mean Lx-Tx relation for the clus- 
ters that we match by means of a linear regression. If we 
look at the positional offset versus offset in L x from the 
mean L x -T x relation at a given T x (given by the quantity 
DL X = Alog(L x /E(z) h^ergs^ 1 )) presented in Fig. 1131 we 
note that the NUV blue BCGs (asterisks surrounded by a 
diamond) tend to cluster at positive values of X-ray excess 
and very small distances from the cluster centres. This was 

10 This is required because bluer galaxies sometimes have an 
overestimated photometric redshift. 
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Figure 12. Positional offset between the BCG and the cluster 
X-ray centre for clusters in common with the ACCEPT sample. 
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Figure 13. X-ray luminosity excess vs. BCG/X-ray positional 
offset for clusters that we have in common with the ACCEPT 
sample. NUV blue (i.e. NUV-r < 5) BCGs are identified with an 
asterisk symbol surrounded by a diamond. There is an obvious 
tendency for star forming BCGs to lie closest to their host clus- 
ter's X-ray peak while the normal red non star forming BCGs are 
the furthest. Moreover, blue BCGs are in clusters that lie above 
the mean Lx-Tx relation. 

the diagnostic used by Bildfell et al. (2008) to infer the link 
between cooling , galaxy position and galaxy blue core. Here 
we are in the position of showing that such a conclusion is 
supported by the analysis of the cooling times (Fig. 1141 up- 
per panel) and entropy (Fig. 1141 lower panel) as provided 
by Cavagnolo et al. (2009). Therefore we can confirm Bild- 
fell et al. (2008) results that a fraction of blue BCGs can 
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Figure 14. Cluster cooling time (upper panel) and entropy (lower 
panel) versus BCG's offset from the the cluster centre. Symbols 
as in the previous figure. NUV blue (i.e. NUV-r < 5) BCGs are 
hosted by cluster with short cooling times and are at low distances 
from the cluster centre. 



be explained by cluster centres where cooling flows supply 
cold gas for star formation. That is, BCGs in low entropy 
clusters may not have enough star formation to be optically 
blue, but they are definitely below the NUV-r red-sequence 
(Fig. E5]). 

In a companion paper (Pipino & Pierpaoli, 2010) we 
further link the presence of a cooling flow (and hence the 
likely presence of a blue BCG) with an enhanced Sunyaev- 
Zeldovich (Sunyaev & Zeldovich, 1970) signal in order to 
study the bias induced in SZ-based cluster catalogues. Fi- 
nally, a more general comparison of our cluster catalogue 
with nearly 1000 clusters observed in the X-rays at a lower 
spatial resolution and available in the literature is presented 



Figure 15. Cluster excess entropy versus BCG's NUV-r colour 
for clusters in common with ACCEPT. 



in the main catalogue paper (Szabo et al., 2010), where the 
correlation between richness and relevant quantities such as 
Lx and Tx is made. 



7 CONCLUSIONS 

In this paper we characterize in terms of magnitude and 
colours a sub-sample of more than 14300 BCGs drawn from 
a parent catalogue of more than 220000 galaxies in 69000 
clusters based on the matched filter method (Kepner et al., 
1999, Dong et al., 2008) applied to the SDSS DR6 (Szabo 
et al., 2010) 

In agreement with previous works, the BCG luminosity 
is found to a have a redshift evolution broadly consistent 
with pure "aging" of the galaxies. Richer clusters tend to 
have brighter BCGs, however less dominant than in poorer 
systems. 4-9% of our BCGs are at least 0.3 mag bluer in 
the g-r colour than the red-sequence at their given redshift. 
Such a fraction does not change if we restrict ourselves to 
clusters with richness above 30, and decreases to the 1-6% 
above a richness of 50. In this case 3% of them are 0.5 mag 
below the red-sequence. A preliminary morphological study 
suggests that the increase in the blue fraction at lower rich- 
nesses has a contribution for the increase in the fraction of 
spiral galaxies. Therefore we suggest a cut at richness of 
100 (and M T < —22.5 mag) if one wants to focus only on 
early-type galaxies. In terms of redshift evolution, the over- 
all blue fraction goes from ~5% in the redshift range 0.1-0.2 
to ~10% in the redshift bin 0.2-0.3. The blue fraction seems 
to increase at higher redshifts, however the scatter in the 
colours and the fact that the catalog is no longer complete 
hamper us from having firm conclusions. We show that a 
colour selection based on the g-r red-sequence or on a cut 
at colour u-r> 2.2 can lead to missing the majority of such 
blue BCGs. Finally, the blue fraction increase by a factor 1.5 
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at most when the study is extended to the three brightest 
galaxies of each cluster. 

The fraction of blue BCGs is in broad agreement with 
previous works (e.g., Crawford et al., 1999, Edwards et al., 
2007, Bildfell et al, 2008) which showed that about one 
quarter of BGCs show emission lines and optical blue cores 
associated with recent star formation; only a smaller fraction 
has the star formation extended (in time and space) enough 
to make their total colour blue. 

We also show two applications of the Szabo et al. BCG 
catalogue. The first extends the colour analysis to the UV 
range by cross-matching our catalogue with publicly avail- 
able data from Galex GR4 and GR5. We show a clear corre- 
lation between offset from the optical red-sequence and the 
amount of UV-excess. Pipino et al. (2009a) showed a one-to- 
one correlation between optical blue cores created by some 
residual star formation and UV-excess in ellipticals. There- 
fore we can infer that 8% of our BCGs are offset from the 
main sequence because of central recent star formation. The 
fraction of BCGs with even lower residual star formation 
that can be noticed only by means of UV-optical colours 
can be as high as 15%. 

We compare the colours of our BCGs in clusters in com- 
mon with the maxBCGs (Koester et al., 2007) catalogue. In 
7% of the cases maxBCG miss the brightest object because 
is more than 0.3 mag away from the red-sequence. In this 
cases, the galaxy classified as BCG by Koester et al. (2007) 
is still one of the brightest (but not our 1st ranked) galaxies 
in our clusters, this means that maxBCG uses the 2nd or 
the 3rd brightest galaxy as seed for their clusters. Such a dif- 
ference is not enough to explain the differences between the 
AMF and the maxBCG catalogues (see Szabo et al., 2010). 

We cross-match our catalogue with the ACCEPT clus- 
ter sample (Cavagnolo et al., 2009), where accurate tempera- 
ture, density and entropy profiles of the Intracluster medium 
can be found. We find that blue BCGs tend to be in clus- 
ters with low entropy, short cooling times. That is, the blue 
light is presumably associated to gas feeding of recent star 
formation by cooling flows (Bildfell et al., 2008). 
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